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The Lindblad approach to open quantum systems is introduced for studying the dynamics of a 
single trapped ion prepared in nonclassical motional states and subjected to continuous measurement 
of its internal population. This results in an inhibition of the dynamics similar to the one occurring 
in the quantum Zeno effect. In particular, modifications to the Jaynes-Cummings collapses and 
revivals arising from an initial coherent state of motion in various regimes of interaction with the 
driving laser are dealt in detail. 
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■ In recent years, several efforts have been made to clarify the role of the measurement process during the evolution 
L_j6f a quantum system. Besides the conceptual viewpoint, this also has a practical relevance in predicting the results 
ly-^of experiments and possibly pointing out some working conditions where genuine quantum phenomena are expected 
Despite of the different languages, ranging from group algebraic and path integral techniques to stochastic 
evolution models |^,|^, which have been employed to analyze the quantum measurement problem, all these approaches 
O^I'ecognize that a measured system is not closed but interacts with the environment schematizing the meter. An 
equivalence loop connecting these different descriptions has been established in |6| , resulting in an unified picture which 
O i ncorporates the measurement process through an unique parameter proportional the strenght of the system-meter 
ifiteraction. In the present paper, we apply this scheme for investigating the influence of a continuous measurement 
(~| process on the dynamics of a trapped and laser-irradiated ion. This system has lately attracted growing interest, 
culminating in the experimental generation of nonclassical motional states j7||^. The possibility to read out the 
^ vibrational state of the ion following the evolution of its internal levels under a Jaynes-Cummings type interaction has 
een at the heart of the experimental procedure. As we shall see, by including the measurement process an inhibition 
^of the two- level transitions dynamics emerges, which could be made observable via the quenching of collapses and 

revivals in a coherent-state prepared ion. 
^ . , Our physical system is a single trapped ion interacting with a classical laser field and undergoing a continuous 
measurement of its internal population. In practice, one is dealing with a collection of independent evolutions of 
" single ions starting from the same initial state and averaged results are actually registered as outcomes of the ex- 
periment. According to the general discussion of we can either model each single history as a realization of a 
pure-state stochastic evolution, then averaging over the associated selective measurement results, or we can consider 
a deterministic ensemble evolution which directly leads to averaged nonselective predictions. Our choice of the latter 
method relies on the exact solvability of the model. The analysis is based on the master equation for the reduced 
density operator p(t) of the system, which results from tracing out the variables of the external environment from the 
overall system -I- environment density operator and is written in the Lindblad form |^,^ : 

^ ^ -^[H{t),m] - [A, Km , (1) 

where H{t) = Hq + Hint{t) describes the dynamics of the unmeasured closed system and the Lindblad operator, 
representing in general the influence of the environment on the system, is related in this case to the measured 
observable A. The parameter k, expresses the strength of the coupling of the measured system to the meter 
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hereafter assumed to be time-independent. The Hamiltonian Hq describes the ID-harmonic oscillator associated to 
the center-of-mass degree of freedom and the two internal states of the ion: 

Ho = H""' + H^^ = fiuj a) a + huj2i , (2) 

being lo the vibrational frequency, UJ21 the two-level transition frequency, d and ir^ the boson annihilation operator 
and the pseudo-spin z-operator respectively. The trap frequency is supposed large enough to neglect the atomic 
spontaneous emission (strong confinement limit H). The interaction Hint{t) with the laser field is modeled by a 
nonlinear multiquantum Jaynes-Cummings model(JCM) which, in rotating-wave approximation, is written as 

= ^ e'^^' cos[77(a + d^) + ^]a-+ H.c. , (3) 

where ujl is the laser frequency, Qq the fundamental Rabi frequency, a- the lowering operator and ry = lul^xsql/c 
the Lamb-Dicke parameter related to the position standard quantum limit of the ion. In Eq. (^) a standing-wave 
laser field is considered, the phase (p fixing the position of the trap potential with respect to the wave. The measured 
observable A depends upon the specific experiment. If the population of the internal ground state is measured, for 
instance by collecting the fluorescence emitted after stimulated transition to a third auxiliary level as done by Meekhof 
et al. 01, the operator A is the occupancy of level |, i.e. A ^ a _a ^ in our formalism. 

It is convenient to introduce the representation defined by the eigenkets of /foj \S,n)^ S'=|,'f,n = 0,...,(X), with 
matrix elements psn,S'm{t) = {S,n\p{t)\S' ,111) . Moreover, let us assume that the ion is in the low excitation regime 
{oj ^ SIq) and the laser is tuned to the A:-th vibrational sideband, i.e. ojl — 0J21 + kuj, k lE Z, two conditions usually 
fulfilled in laboratory [Q. This implies that transitions between states | i,n) and | t,rt -I- fc), involving the exchange 
of k vibrational quanta, are resonantly enhanced, whereas all the off-resonant couplings are rapidly oscillating with 
frequency lu and can be disregarded (JCM approximation |^). As a notable case, when if = ±7r/2 and 77 <C 1 
(Lamb-Dicke limit), the well-known linear one-quantum JCM or anti-JCM operators ri{aa+ + d^(7_), ri{d''a+ + da-) 
are recovered by tuning to the first red or blue sideband, k = —1 01 k ~ 1 respectively [ p^ . In what follows we focus 
on the properties of the internal dynamics of the ion, which can be extracted from the reduced atomic density matrix, 
o'SS'it) = J2'^=o PSn-^S'nit). We can therefore restrict ourselves to the subset of the density matrix psn,s'm having 
n — m. According to the master equation (|l|) and for a generic fc-quantum resonance, their time development is ruled 
by: 

PinTn+fe(i) = (it^21 — PlnTn+k{t) — i^nn+kf^ [Pln+kTn+k{t) — Plnln{t)] 6*""^ 



Plnln{t) = -iflnn+k/^ [p-f„+fex„ (i)e*"-f'* - px„|„+fe (i)e" 

PTn+k^n+k{t) = —iflnn+kf^ [Plntn+k{t)e~'^'^^* — P]n+k. 
Pin\n+k{t) = {iuJ2l — k/2) Pin^n+k{t) — i^nn+kf^ [Pin-i 
P]n+klnit) = { — iuj21 ~ k/2) P]n+kln{t) — i^nn+k/'^- [Pinln{t) — P\n+k\n+k(^)\ C"*""^* 



where, as usual, the nonlinear /c-quantum Rabi frequencies expressed in terms of the Laguerre polinomials C}^ 

f I \ 1/2 

f^„„+fc = f^o(n| cos[77(a + d^) ^]|n + fc) = fio [c'^ + (-l)'=e-''^] {i^f [j^^^^ e'^" Ci{ri^) (5) 



have been introduced ||9[. In the Lamb-Dicke regime (ry ^ 1) they reduce to the linear ones, e.g. f2„„+i — •qO.Q^n -\- 1 
for one-quantum interaction. The following solution of (^) holds for n > — fc: 

Plnln{t) = ]^{j)]n+k]n+k + Pini«) (0) - {p]n+k^n+k " PirUn) (0) COs(u;„„+fci) 
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Sm{Wnn+kt) + tlm(pxn|„+fc(0)) sm{wn„+kt) 



'^nn+k 



P]n+k]n+k{t) — (^P]n+k]'n+k + Pining (0) - Plnln{t) , (6) 

PinTn+fc(0 = e^'''"*~'^*/4(Re(pi„T«+fe(0))e""*/^ + lm(pi„t„+fe(0))[cos(u;„„+fci) 



Sin('u;„„+fet) - i(pi^n+k\n+k - Plnln) (Q) sin(-Wnn+fc^) 



where 
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Wnn+k — \j ^nn+k 



(7) 



and pin+kin{t) = Px„|n+fc(^)- (H)' holding within every fixed Jaynes-Cummings manifold | i,n), | 1,71 + fc), are 

formally equivalent to Eqs. (84)-(85) of Rcf. describing a two-level system subjected to a continuous occupancy 
measurement. When n — the effect of the measurement disappears and oscillations with angular frequency flnn+k 
result for fixed n. In the opposite limit of strong measurement coupling the frequency becomes imaginary and 
an overdamped regime occurs in which transitions are inhibited, the so-called quantum Zeno effect. The borderline 
between the two regimes is ruled by the threshold Wnn+k = 0, corresponding to a critical measurement coupling 
^rin+k ~ 4f^nn+fc- In & general situation where different vibrational subspaces are involved, the n-dependence of this 
critical value globally results in a smoother transition. The measured occupancy Pi(t) — cult) is straightforwardly 
evaluated from Eqs. (^) once the initial conditions are specified. We assume here, according to |^], that no entangle- 
ment between the internal and motional degrees of freedom occurs and only the ground internal level is populated, 
yielding 
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the matrix elements characterizing the number state distribution of the vibrational motion. The presence of 
entanglement, crucial to extend the description to Schrodinger cat states can be included as well. Eqs. (||) 
displays the effect of the measurement for an arbitrary vibronic distribution, namely an overall damping of the signal 
amplitude and, for each vibrational component, both a frequency shift ^ and a k- weighted term. 

Let us now specialize the internal state evolution (g) to an initial coherent distribution, /9^™(0) = n"e~"/n!. In 
the limit in which no measurement is present (k = 0) and a linear Jaynes-Cummings interaction is realized by a 
sufficiently small Lamb-Dicke parameter, Pi [t) undergoes collapses and revivals, a well-known quantum phenomenon 
due to dephasing and rephasing between the various Rabi oscillations, as firstly predicted by Eberly et al. and 
recently observed in cavity QED by using Rydberg atoms Fig. la shows the collapses and revivals occurring 

in a quantized trap, the center-of-mass being in a coherent state of average quantum number n. The modifications 
intervening when the same system is coupled to the meter are depicted in Fig. lb for a strength k = 10~^k3"*, in 



units of a reference value 
(1) leads to 
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Provided k <C k-'^'* 
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Int(n), an approximate evaluation of the series in 



Pl{t) ~ 1 + cosiijnoVn +lt) exp ( - '-t^ 
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which holds for a timescale t <C 2^/^/rjQ,Q, still including all the collapse evolution, and generalizes the result already 
quoted in [ p| for k = 0. Thus, even in this weak coupling regime, the measurement has a leading influence at short 
times due to its linear dependence in the exponent of the envelope in (^) . Similar considerations can be repeated for the 
nonlinear multiquantum Jaynes-Cummings interaction, which is restored outside the Lamb-Dicke limit and has been 
extensively studied by Vogel and de Matos-Filho in case k = Two representative examples, corresponding to 
one-quantum and two-quantum couplings are shown in Figs. 2 and 3 respectively The parameters have been chosen to 
reproduce, in the closed-system limit, the experimental results reported by |Q for a nonlinear one-quantum interaction 
and the same reference ratio k/kq^** = 10~^, fc = 1, 2, has been maintained as before. In all the cases the measurement 
tends to wash out collapses and revivals making more difficult to infer the quantized nature of the vibrational motion 
in the trap. For larger values of /i/^p"* the suppression is even more pronounced, while for > 1 a complete 

freezing of the population dynamics to the initial occupation is achieved, regardless of the initial vibrational state. 
A qualitatively similar quenching of the collapse-revival phenomena for the linear one-quantum Jaynes-Cummings 
dynamics of a single two level atom in an optical cavity is described in by including both spontaneous emission 
and cavity damping mechanisms. The relative insensitivity to spontaneous emission damping reported there further 
enforces the possibility of neglecting such effect in our description. 

Trapped ions seem particularly adequate to investigate the vanishing of revivals in coherent states induced by 
the measurement of their occupancy probability. Relaxation mechanisms due to spontaneous decay of the electronic 
transition or mechanical dissipation of the ion vibrational motion are indeed negligible |jl^,|l^ . A technical source of 
decoherence has been instead identified in Q and ascribed to laser and frequency trap fluctuations, contributing to 
a finite Hnewidth for the electronic transition depending upon the vibrational state, 7„ — 7o(ri + 1)" '', 70 — 11.94 
kHz. This yields an average decay time = 7^^ = 31.2 /is for n = 3.1. On the other hand the data collected in 



the quantum Zeno experiment ||lj 
are fitted through Eq. (m) with k 



In 

exploiting the same auxiliary transition ^ ^P?,/2 of ^Be"*" 

4.9 • lO^'s^^ [1). The corresponding decay time is t = 4k^^ 



also used in [Q, 
816 /iS, one order 
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of magnitude larger than the technical decoherence decay constant, suggesting that an improvement in the stability 
performances of the experimental setup by roughly a factor 10^ should imply a measurement-dominated decoherence. 
The previous numerical example corresponds to a ratio k/ k'^^^j^^. = 1.1 • 10~^, i.e. a weak coupling regime. An increase 
of this ratio by two orders of magnitude, for instance attainable by using the same measurement configuration but a 
smaller Rabi frequency fio, allows one to test the region where neither a closed system evolution nor a von Neumann 
instantaneous projection of the state are suitable for reproducing the experimental results. A description in terms of 
the Lindblad equation (|l|) is instead mandatory. 

The predicted effect does not exhaust the influence of the measurement process on the ion dynamics. As inves- 
tigated in in a different context, the measurement also indirectly affects the average vibrational motion of the 
trapped ion, a quantum damping without classical analogue whose features will be described in [^ . Furthermore, 
analogous decoherence effects can be studied within the same framework in other physical systems, atomic Bose- 
Einstein condensates being particularly appealing in connection with a recent analysis of their collapses and revivals 
dynamics [Q. From a general viewpoint, experimental investigations of the effect of the measurement process as 
the one proposed here will allow for a quantitative description of this ultimate source of decoherence, a topic which 
besides its established importance in defining the boundary between quantum and classical worlds p9[ | is becoming 
crucial for the issue of quantum computation pU. 
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FIG. 1. (a) Ground state occupation versus time for a linear one-quantum resonance and an initial coherent state distribu- 
tion, showing collapses and revivals. The parameters n = 3.1, 77 — 10~^ and f2oi/27r = 4.75 KHz have been used and k = 0. 
(b) Same as in (a) except k/kq"* — 10~^. 
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FIG. 2. (a) The same as in Fig. la but for an one-quantum nonlinear Jaynes-Cummings model, Q,q\/2'k = 94 KHz and 
77 = 0.202. (b) K/«§r* = 10-^ 

FIG. 3. (a) The same as in Fig. 2a but for a two-quantum nonlinear interaction, f2o2/27r = 13.4 KHz. (b) k/kq"* = 10~^. 
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